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Abstract 

We consider the role of the velocity in Lorentz-violating fermionic quantum theory, 
especially emphasizing the nonrelativistic regime. Information about the velocity will be 
important for the kinematical analysis of scattering and other problems. Working within 
the minimal standard model extension, we derive new expressions for the velocity. We find 
that generic momentum and spin eigenstates may not have well-defined velocities. We 
also demonstrate how several different techniques may be used to shed light on different 
aspects of the problem. A relativistic operator analysis allows us to study the behavior 
of the Lorentz-violating Zitterbewegung. Alternatively, by studying the time evolution of 
Gaussian wave packets, we find that there are Lorentz-violating modifications to the wave 
packet spreading and the spin structure of the wave function. 
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1 Introduction 



Recent work has stimulated a great deal of interest in the possibility of there existing 
small Lorentz- and CPT-violating corrections to the standard model. If violations of these 
fundamental symmetries exist at low and medium energies, they could represent powerful 
clues as to the nature of Planck scale physics. Within the context of effective field theory, 
a general Lorentz- violating extension of the standard model has been developed ^1 12] , and 
the stability and renormalizability |3| of this extension have been carefully studied. 
However, the general standard model extension (SME) is extremely complicated, and 
even superficially simple questions about its physics may have subtle and even ambiguous 
answers. For example, the study of the gauge invariance properties of and finite radiative 
corrections to this Lorentz- violating field theory has proven to be a fruitful source for new 
theoretical insights El El El • Many other elementary questions about the SME still 
remain unanswered. 

The SME provides a framework within which to analyze the results of experiments test- 
ing Lorentz violation. To date, such experimental tests have included studies of matter- 
antimatter asymmetries for trapped charged particles (TUl HH 1121 US] and bound state 
systems [T^ [T3] . determinations of muon properties ^Hl analyses of the behavior 
of spin-polarized matter p!Hl ITU' , frequency standard comparisons 1201 1211 122] , measure- 
ments of neutral meson oscillations |i2Sl 1211 12^1 , polarization measurements on the light 
from distant galaxies [201 123 l2Hj , and others. The analysis of the relevant experimental 
data requires a good understanding of the behavior of elementary particles in the pres- 
ence of Lorentz violation. However, there still remain many aspects of Lorentz-violating 
quantum mechanics about which much more could be known. 

In this paper, we shall examine one important aspect of Lorentz-violating fermionic 
quantum mechanics: the behavior of the velocity. A detailed understanding of the role 
of the velocity in Lorentz-violating physics is important for analyses of Lorentz-violating 
scattering processes. The relevant kinematics of such processes may depend sensitively on 
the character of the velocity. For example, the question of whether the vacuum Cerenkov 
reaction [2211201 e~ —>■ e~7 can occur is obviously intimately related to whether or not the 
initial electron's velocity is superluminal. There have been some prior investigations into 
the properties of the velocity in the presence of Lorentz violation. Some of these analyses 
have studied the effects of specific forms of Lorentz violation, involving selected terms from 
the SME or Lorentz violation through "double Special Relativity" [121E3ISlj- There 
has also been some analysis [35j of velocities in nonrelativistic Lorentz-violating theories 
(obtained from their relativistic counterparts through the use of the Foldy-Wouthuysen 
transformation [2^1) • However, there has been no comprehensive treatment of the velocity 
in the context of the Lorentz-violating SME. 

This study is a natural outgrowth of earlier work on the effects of Lorentz violation 
on scattering. Previously, we have looked at the impact of specific Lorentz-violating 
coefficients on particular scattering processes [2ZlEHj- Lorentz violation affects particles' 
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velocities and hence the kinematics of reactions. These effects can be just as important 
as any change in the dynamics. Exphcit relations between the velocity, momentum, and 
spin are required for these types of calculations. To date, a systematic treatment of these 
relations is lacking. It is the purpose of this paper to fill this gap and to obtain a more 
detailed understanding of the role of the velocity in Lorentz-violating quantum theories. 

In this paper, we shall concentrate primarily on the nonrelativistic regime. There are 
several reasons for doing this. First, this regime has thus far received little attention; 
previous kinematical studies have tended to focus on relativistic speeds and particular 
systems, such as ultra-high-energy cosmic rays. Second, only nonrelativistic quantum 
theory possesses a satisfactory probabilistic interpretation, and we shall need to make 
use of the probabilistic method when we analyze the time development of wave packets. 
Third, consideration of the nonrelativistic case will help us to distinguish which of the 
theory's complications are due to the effects of Lorentz violation and which are merely 
byproducts of the matrix structure of Dirac theory. For example, the velocities we en- 
counter will generally be spin-dependent, but in the relativistic matrix theory, only the 
helicity component of the spin is a conserved quantity. If we restrict our attention to the 
nonrelativistic limit, in which all the spin projections are time- independent, we may avoid 
this difficulty. 

Although we shall mostly be interested in particles with nonrelativistic speeds, we 
shall use a number of inherently relativistic methods in our study of the velocity. The 
velocity operator in a relativistic fermion theory has a complicated matrix structure, 
and we shall examine how this structure is modified by the presence of Lorentz-violating 
parameters. This analysis will allow us to address questions about aspects of the theory 
that are inherently relativistic, such as Zitterbewegung. We shall also use fully relativistic 
energy-momentum relations to calculate group velocities. Finally, when it is possible for 
us to display concise exact expressions for the velocity, these expressions will necessarily 
be relativistic. 

In keeping with our nonrelativistic viewpoint, most of our considerations will apply 
only to the case of massive particles. Moreover, we shall generally consider only one 
Lorentz-violating interaction at a time. This is a sensible approach if all the Lorentz- 
violating parameters are small, so that we need only work to leading order in each param- 
eter. In the presence of multiple forms of Lorentz violation, the various Lorentz-violating 
contributions to the velocity are then simply additive. However, although we shall concen- 
trate on obtaining the leading order corrections to the velocity, we shall, as noted above, 
also present exact results when possible. An exact treatment offers substantial additional 
verification that the effects we uncover are indeed meaningful. 

We shall introduce the Lorentz-violating coefficients relevant to a purely fermionic 
theory in Sectional Then we shall immediately specialize to a theory containing only one 
specific Lorentz-violating parameter: a timelike axial vector term. Such a term generates 
many interesting effects, most of which may be analyzed exactly. We shall begin our 
study of the velocity with an analysis of the propagating modes of the fermion field in 
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the presence of this term. From the energy-momentum relation for these modes, we can 
extract a group velocity. We then move on to a general analysis of all the possible Lorentz- 
violating terms, using the relativistic Dirac algebra (Sectional). We shall determine the 
Zitterbewegung-free contribution to the velocity operator, which governs the bulk motion 
of the particles. Finally, we turn in Section 0] to an examination of the time evolution of 
Gaussian wave packets. This wave packet analysis will allow us to resolve several subtle 
questions about the role of the velocity that will arise in the course of our discussions. 
Finally, we shall summarize our conclusions in Section |S1 

2 Lorentz- Violating Fermion Theory 

2.1 Lagrangian Structure 

Since our studies of Lorentz violation will focus only on kinematics, we shall work with 
a theory that describes a single species of noninteracting fermions. Including only the 
minimal, superficially renormalizable Lorentz-violating terms, the Lagrange density for 
this theory is 

C = iji^T'^d^ - M)^, (1) 

where 

M = m+^- + -H'^'a^, + im^'j^, (2) 

and 

r'^ = 7^ + c-"-!. - d-^-J^J, + e'^ + + ^g^^'^a,,. (3) 

These represent the only superficially renormalizable couplings that are possible in a 
purely fermionic theory. However, some of the couplings are more interesting than others. 
In particular, ms does not violate Lorentz symmetry, and it may be absorbed into the other 
coefficients by means of a field redefinition [3^] . We shall assume that such a redefinition 
has already been performed and set = 0. Moreover, e, /, and g are inconsistent with 
the coupling of the fermion field to standard model gauge fields. However, we shall include 
these three coefficients for the sake of greater generality. 

2.2 Propagation Modes in the Presence of b'^ 

We shall begin our analysis by considering the specific example of a 6-type interaction. 
This form of interaction is particularly illustrative, and, moreover, the results of our 
careful analysis of this theory will be needed when we consider the more complicated 
effects of d'^'^ in the Appendix. The quantization of this particular system has previously 
been examined in j3U], where some properties of the velocity were noted. In general, 
the spacetime direction of b is arbitrary, but if b is timelike, we may consider it in an 
observer frame in which its three spatial components 6* vanish: 6^ = {W,0). Since the 
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time dimension plays a special role in the canonical quantization of the Dirac field, this 
is the most natural frame in which to quantize the theory. 

In order to quantize the 6-modified theory, we must determine the free propagation 
modes of the fermion field. So we shall solve the free momentum-space Dirac equation, 
with the effects of the b term included. Because the matrix 75 features prominently in 
the theory, we shall use the Weyl chiral representation for the Dirac matrices: 
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We consider a fermion mode with energy E and three-momentum p 
by 7°, the Dirac equation for this mode may be reduced to 



E - 6° + P3(T^ 



—m 



—m 



E + b^ - psa^ 



u{p) = 0. 



(4) 



p^z. By multiplying 



(5) 



If the fermion's spin is quantized along the 2;-axis, with eigenvalues ±|, then we may 
replace —>■ s. The eigenvalue condition for E then becomes 



E' 



m 



+ (sps - b^y = rr? + {s\p \ - 6°)' 



(6) 



Note that even though there is no breaking of rotation invariance, the energy depends 
upon the spin direction, through the helicity s. The explicit electron spinors corresponding 
to a specified momentum are given in [SHj; the corresponding positron spinors are precisely 
analogous. These solutions may easily be generalized to describe modes with arbitrary 
three-momentum, so long as the spin is quantized along the direction of the motion, so 
that s represents the helicity. In the nonrelativistic limit, both the positive- and negative- 
energy spinors approach their usual limits, provided that 6° is small compared to the 
mass. 

A number of odd facts follow from the dispersion relation (jH)). These properties have 
been previously noted in [3 EI], but we reiterate them here, because a proper understand- 
ing of these somewhat counterintuitive results is necessary if one is to possess a complete 
understanding of the meaning of the velocity in Lorentz- violating physics. Restricting our 
attention to the E > Q modes, we see that the p = Q mode has energy ^Jm? + (6°)^, but 
that this is not the lowest-possible energy. Indeed, for s|p | = 6°, the energy is only m. So 
for 6° 7^ 0, a particle with p = can release energy by absorbing momentum. 

However, such a particle, with vanishing three-momentum, is not actually stationary. 
Because of the Lorentz-violation, the group velocity for a wave packet centered around 
p = is nonvanishing. In fact, the group velocity for an arbitrary wave packet that is 
well-localized in momentum space is 



V^E 



p 



sb^p 



E 



(7) 
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Note that both the numerator and denominator contain contributions that are first order 
in }p. For vanishing p and IP < 0, the velocity points along the direction of the spin. 
Since Vg is a vector and the spin is an axial vector, this result represents a clear physical 
manifestation of the parity violation arising from 75. Moreover, if the spin is not quantized 
along the direction of p, then the wave function will contain a superposition of two different 
energy eigenstates, with the two energies corresponding to different group velocities. (It 
is perhaps somewhat unsurprising that there are momentum eigenstates which do not 
possess a unique group velocity, because the very concepts that we use to define the 
velocity — the flows of energy and particle density — lose some of their usual meanings 
when breakdowns of Lorentz symmetry occur.) This also creates the possibility that the 
wave packet may bifurcate, separating itself into two distinct spin states. However, since 
the two spin components perpendicular to the momentum are not constants of the motion 
in the Dirac theory, neither is it clear that such bifurcation will in fact occur. Only a 
careful analysis of the problem in terms of wave packets can answer this question. We 
shall perform such an analysis in Section 14.21 for a wave packet with vanishing mean 
three- momentum, we find that, when Ip is small compared with the momentum spread 
of the wave packet, no bifurcation of the wave packet in the plane perpendicular to the 
initial spin can occur. 

3 Operator Methods 

Because the group velocity for a wave packet may, under certain circumstances (as in 
systems displaying anomalous dispersion) become meaningless or even undefined, some 
further confirmation that ((Tj) represents a real physical velocity is desirable. This confir- 
mation may be found through an analysis of the operator structure of the Dirac theory. 

In Dirac theory, the velocity possesses a number of decidedly nonclassical properties. 
In the presence of Lorentz violation, this fact is even more evident. Indeed, the Lorentz- 
violating contributions to the velocity may depend upon operators that have no classical 
analogs. Moreover, unlike ordinary oscillatory Zitterbewegung terms, the Lorentz- violating 
additions to a fermion's velocity will unavoidably affect the bulk motion of that particle. 

We shall examine the effects of all the Lorentz-violating interactions included in 
in this way, beginning with the M terms and, in particular, the specific M = m — 6*^7*^75 
considered in Section 12.21 After obtaining an operator result analogous to ((Tj) for the 6° 
theory, we shall generalize our methods to deal with the remaining M terms and the F'^ 
terms as well. 

3.1 Potential (M) Terms 

When M = m — b^'y^'j5 and F^ =7^, the Heisenberg equation of motion for the position 
operator Xk is 

Xk = i[H,Xk], (8) 
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where 

H = ajPj + 13m- 756° (9) 

is the single-particle Dirac Hamiltonian. The Dirac matrices are, as usual, aj = 'y^'j^ and 
(3 = 7°. So the time derivative of Xk is simply ak, just as in the Lorentz-invariant theory, 
and a remains the correct operator to represent the velocity. 

We shall therefore consider the equation governing the time evolution of at- This is 

dfc = t[H, ak] = I {-2akH + 2pk - 2lPaki^) ■ (10) 

When 6° = this has the exact solution (keeping in mind that p and H are constants of 
the motion) 



ak{t)=PkH ^+ ak{0)-pkH ^ 



(Ill 



When 6° 7^ 0, finding the solution becomes much harder. However, we can determine 
the solution to first order in W if we consider only the case in which p = 0; this is again 
permissible because p is conserved. (Note that we are not boosting the theory into the 
p = frame, because this would generate a spacelike component for b. Instead, we are 
merely restricting our attention to the p = subspace of the theory.) In this case, the 
time evolution of a is given by 



akit) ^ + 

m 



6° 

«fc(0) + -T.kl3 
m 



e-^'"K (12) 



That (fT^ is indeed correct to first order in Ip becomes quite clear when we observe several 
facts about the p = subspace. The spin operator S = 075 is a constant of the motion, 
and /? is (9 [(6*^)^], when these operators are considered within this subspace. Moreover, 
when p = 0, H = (3m + (9[(6'')^], so the time evolution of the 6-linear terms is governed 
entirely by their commutators with (3. 

The second term on the right-hand side of (fT^ contains the usual Zitterbewegung os- 



cillations with frequency u = 2\H\ = 2^m'^ + (W)'^. However, the Lorentz violation may 

modify the structure of this Zitterbewegung. The eigenvalues of the operator ak{0) + ^T,k/3 
(as well as those of H) essentially determine the length scale over which the Zitterbewe- 
gung occurs. The eigenvalues of afc(O) + ^S^/? are iy'l + {W/my, and since these have 
no 0{b^) contribution, we see that there is no leading-order modification of the length 
scale over which the Zitterbewegung motion occurs, at least for particles with vanishing p. 
We do see indications of an effect — a lengthening of the scale — at higher order. However, 
this effect would be exactly compensated by the increase in the frequency of the oscilla- 



tions from 2m to 2^m'^ + (6°)^. In any case, since we have not conducted a systematic 
expansion to (9[(6°)^], we cannot make a definitive statement about the higher-order cor- 
rections. If there were a correction to the length scale of the Zitterbewegung, the effect 
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would manifest itself, for example, in the Darwin term in the Hamiltonian. The Dar- 
win term is generated by Zitterbewegung and is proportional to the Zitterbewegung scale 
squared. 

However, the first term on the right-hand side of (jl2j) is more interesting. It has 
precisely the form that we would expect on the basis of ((Tj). It does not oscillate and 
is not related to interference between positive and negative energy solutions, so it will 
not be swamped by the ordinary Zitterbewegung, nor can it be eliminated by considering 
wave packets with only positive-energy components. On the basis of this result, we may 
conclude that the velocity Vg is not merely an artifact of the approximations used in the 
calculation of the group velocity; it is a real effect deriving from the peculiar relativistic 
quantum mechanics of this system. 

The result ()12|) also demonstrates that velocity and momentum are fundamentally 
different objects in this theory. It is already familiar from the ordinary Dirac theory 
that the velocity operator has a number of unattractive features, and ai do not 
commute with one-another for k ^ I, nor do they commute with the Hamiltonian, and 
the only eigenvalues of are dzl. However, if we choose to restrict our attention to wave 
packets containing only positive energy components or to ignore the rapidly oscillating 
part of ak{t) (which averages to zero), we can sensibly identify the k-th component of 
the velocity with the first term in (fTT|) . pkH~^, which is exactly the classical velocity 
for a relativistic particle. However, in (|12j). there is a new, spin-dependent velocity term 
that exists even at zero momentum. This term can have no classical analogue, because 
of its nontrivial matrix structure. There is therefore no way that we can extend the 
classical correspondence between momentum and velocity to this case. Moreover, it is 
the unavoidable matrix structure of the velocity that leads to the existence of momentum 
eigenstates that are not eigenstates of V . 

We may adapt the approximate solution (jl2|) for ak{t) to account for a much more 
general M. However, we shall see that, of the possible terms included in M, only a\ 

and H^^ generate nontrivial changes to a at leading order. (Beyond leading order, 
any Lorentz-violating interaction that modifies the energy-momentum relation will cause 
a corresponding change in the classical velocity term pkH~^.) In order to derive our 
generalization of (|12|). let us consider a single-particle Hamiltonian whose p = restriction 
is if = jSm + H' . H' is the source of the Lorentz violation, which is parameterized by 
some small quantity e: H' = 0{e). First, let us note that if H' commutes with both ak 
and (5, then H' is a constant of the motion that does not contribute to ak- Of the sixteen 
possible Dirac matrices, four — /, ejki<y^\ and 7-^75 for j 7^ k — have this property. So if H' 
is any of these, the p = solution for ak is 

ak{t) = afc(0)e-2*'3"^*, (13) 

and this result is exact. 

Otherwise, if either [H', j3] or [H', ak] is nonzero, the equation of motion for ak becomes 

Ok = —2iakH + i{H, at} = —2iakH + i{H', ak}- (14) 
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We consider a possible solution to (fT^ . valid through (9(e) and similar in form to (fT^ : 



«fc(t)~V^fc + K(0)-rfc]e-2^^*, (15) 

where is C(e), and 14 is C'(e^)- Inserting this expression into (fT^. we see that the 
solution is valid to (9(e) if 2iVkH = i{H', ak}. Therefore, up to corrections of (9(e^), 

Vk^l {H', «4 ^ {H', ak} p. (16) 
2 Im 

We must now check the consistency of the statement that = (9(e^). This amounts 
to checking that [{H',ak}, P] = 0. We have already considered all the Dirac matrices 
that commute with both and /?. Therefore, since [{H',ak}, P] = —[{//',/?}, a^], and 
any two of the basic Dirac matrices either commute or anticommute, all twelve of the 
remaining matrices — /?, 7-', aj, a^^, Pj^, 7^75, and 75 — satisfy the condition. However, 
only those of them that commute with will contribute to 14 . So it is precisely those 
matrices which anticommute with P and commute with that will ultimately affect the 
velocity. 

In fact, in terms of a general H' = P{M — m), the leading-order contribution to the 
velocity is 

}fi f.k TTOj 

V, = --E,P--P+t^e^'%. (17) 
m mm 

The presence of P in the a and b terms has a simple explanation. To this order, P is 
simply the sign of the energy, and so it is the operator that determines the direction of 
the motion. 

The three terms in (jl7p exhaust essentially all the possible contributions to the velocity 
that can be composed from three- vectors only. At p =0, the only three- vectors that can 
appear in V are the spin, an externally prescribed vector, and a cross product of the spin 
and a prescribed vector. Each of these types is represented in (fT7|) . We also note that, 
like the 6° contribution, the and H'^^ terms in Vk are exactly what could be expected 
from the energy-momentum relations. In the presence of a only, the energy becomes 

= m'^ + {p - df, (18) 

and with only an H^^ term, the energy is 



E^ = m' + p^ + ± 2i7°Vpi+pi, (19) 

where the sign of the square root depends upon the spin orientation. Each of these gives 
rise to a group velocity in agreement with the corresponding term in ()17p . 

If we repeat our earlier analysis of the Zitterbewegung for a'^ and H^^ , we find that 
neither of these makes a first-order contribution to the Zitterbewegung scale. However, the 
higher order cancellation we observed does not occur in these cases. The eigenvalues of the 
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Zitterbewegung component of the velocity depend on the direction under consideration, but 
the Hamiltonian does not, so the motion could take on an elliptic, rather than spherically 
symmetric, character. (Again, however, a higher-order calculation would be required in 
order to verify this.) 

We must now discuss the generalization of (jl7|) to the p case. There can exist 
additional complications in the interpretation of the velocity in this situation; for example, 
as we saw in the theory containing the W term only, there is no unique group velocity 
unless the spin and the momentum are collinear. However, to the extent that we may 
approximate {H, a^} to be a constant of the motion, we may also approximate 



ak{0)-^{H,ak}H-' 



e-^'"K (20) 



For the iP-orAy theory, ()2()|1 is exact for the component of a that points parallel to p, but 
it need not even be a good approximation for the two perpendicular components. 

Fortunately, however, a complete generalization of the operators methods leading to 
()17|) is unnecessary. We may instead utilize observer Lorentz symmetry to determine the 
velocity when the three-momentum is nonzero. Since we have considered the most general 
possible M, we may always boost a particle under consideration into a frame in which 
p = (although it would now be a misnomer to refer to this as the particle's "rest frame") 
and then use p7|) to determine V in this frame. An inverse boost then gives the velocity 
in the original frame. This boosting procedure does not require prior knowledge of the 
particle's velocity, because we are boosting into a frame in which p = 0, rather than into 
one where the particle is stationary. (However, since the velocities involved are ±p/E, a 
knowledge of the Lorentz- violating expression for the energy is required.) This therefore 
concludes our analysis of the M terms. 



3.2 Kinetic (F) Terms 

We now turn our attention to the Lorentz- violating terms present in F'^. The canonical 
quantization of the fermion field requires some care when the c, d, e, f , or g coefficients 
are nonvanishing, because C may include additional time derivative terms beyond the 
usual one. A matrix transformation ijj — >■ Atp may be required in order to ensure that 
po _ ^0 ^j^^ explicit expression for the required A, valid to all orders in the Lorentz 
violation, is given in [H].) In this section, we shall assume for simplicity that any such 
necessary transformation has already been performed. This amounts to setting d^^ = 
^uo = gO _ yo _ g\M _ g However, once such a transformation has been made, we 
may not freely boost the theory into another frame; such a boost would reintroduce the 
troublesome contributions to F''. We will therefore be restricted to considering any theory 
with a nonstandard F'^ in only a single observer frame. 

In this section, we shall only consider the terms c, e, and /. The d and g terms, 
because of their dependences on 75, lead to more complicated calculations; the velocities 
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will necessarily be spin dependent, and this leads to awkward initial conditions. However, 
for the theoretically better motivated d term, a complete analysis is given in the Appendix. 
We shall not consider any Zitterbewegung effects in this section, because simple eigenvalue 
estimates of the distance scales involved do not work when the momentum p is nonzero; 
even in the absence of Lorentz violation, it is very difficult to estimate the scale of the 
Zitterbewegung motion directly for a fermion with nonzero momentum. 

In fact, the velocity in the presence of any of c, e, or / may be found exactly. We 
begin with the fairly straightforward c terms. In the presence of c'^^, the single-particle 
Hamiltonian becomes 

H = ajPj - cijaipj - CojPj + f3m, (21) 
so the time derivative of the position operator is 

ik = oik- cikOi - cofc. (22) 

The velocity and a are no longer one and the same, but the two remain very closely 
related. 

We may solve for ak with the same methods as we used previously. The equation of 
motion 

dfc = i [-2ak{H + cojPj) + 2pk - 2ckjPj] (23) 

has the exact solution 

ak{t) = {Pk - CkjPj) {H + CojPjY^ + 

(24) 

The Zitterbewegung-iree contribution to the velocity is therefore 

Vfc = {pk - CkjPj - CjkPj + CjkCjipi) {H + CQjPj)'^ - cok- (25) 

Note the presence of the inverse Hamiltonian. Whenever H appears in Vk, we must, as we 
noted following equation ((7j), account for the Lorentz- violating modifications of the energy 
eigenvalues when determining the leading-order corrections to the velocity. To first order, 
only the symmetric part of Ckj contributes to V; the antisymmetric part corresponds at 
this order merely to a change in the representation of the Dirac matrices, which should 
have no physical consequences. 

In the presence of or interactions, with Hamiltonians 

H = ajPj + [3{m — CjPj) (26) 

or 

H = ajPj + pm- ip-f^fjPj, (27) 

the solution ()11|) for ak{t) remains unchanged. The only Lorentz- violating contributions 
to V come from the modified relations Xk = ak — (3ek and Xk = ak ~ iPj^fk- The 



"fc(0) - {pk - CkjPj) {H + CQjPj] 



-2i{H+cojPj)t 
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Zitterbewegung-free contributions to Pit) and (/375)(t) in the presence of e and /, respec- 
tively, are (m — ejPj)H~^ and ifjPjH^^. So the Zitterbewegung-free velocities are 

Vk = PkH-^ - Ckim - ejPj)H-^ (28) 

and 

Vk=PkH-' + fkif,Pj)H~\ (29) 

Note that the expression (|29|) for the velocity in the presence of / contains no first- 
order contributions. These exact results are again entirely in keeping with the energy- 
momentum relations E'^ = [m — e ■ pY + p"^ and E"^ = + + (/ ■ p)^- 
This concludes our analysis of the Dirac algebra relating to the velocity. 



4 Wave Packet Analyses 

The study of the Dirac algebra associated with the velocity can be quite illuminating, 
but there remain some questions that can be more satisfactorily answered through the 
use of other techniques. In particular, we have not dealt with the issue of wave packet 
bifurcation that was raised earlier. There are also problems relating to wave packet 
spreading, such as whether ordinary spreading might swamp the Lorentz- violating effects 
that we are considering. It will also be instructive to examine how the spreading itself may 
be modified by the Lorentz-violating parameters. To study these aspects of the problem, 
we must construct and study the particle wave packets directly. 

For our wave packet analyses, we shall use the nonrelativistic Foldy-Wouthuysen for- 
mulation of the theory [3^1 02] • We do this, in part, because only the nonrelativis- 
tic quantum theory has a completely consistent probabilistic interpretation. Also, in a 
nonrelativistic theory, with m large compared to any energy scale associated with the 
Lorentz- violation, we do not need to include the effects of the negative energy modes; this 
eliminates the troublesome Zitterbewegung, whose structure we have already considered 
using the operator formalism. Moreover, since the positive-energy spinor M*(p) is given 
approximately by 

rnxs 



u'{p) 



mxs 



(30) 



in the nonrelativistic limit, we shall only need to consider two-component spinor wave 
functions. 

In Section 14. ![ we shall derive an expression for the velocity operator in the non- 
relativistic Foldy-Wouthuysen representation. We shall then examine how the modified 
velocity affects the bulk motion and spreading of an initially Gaussian wave packet. How- 
ever, we shall ignore all spin-dependent effects, except those which contribute directly to 
the mean velocity. Subsequently, in Section we shall treat the problems related to the 
spin more carefully, examining the possibilities for spin flips and spin-driven wave packet 
bifurcation. 
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4.1 Wave Packet Velocities and Spreading 

We begin with the nonrelativistic (Foldy-Wouthuysen) single-particle Hamiltonian, as 
derived in j32]-^ This is 



H 



FW 



2m 



+ 



PjPk 



+ -bj + mdjo - -mejM9kio + -(^juHki + [aj - m (cqj + 90) - me^ 



Pj 



m 



bo6jk — m {dkj + dooSjk) — rntkim y^^mij + gmoo^jij — ^jkiHio 

+ I m {doj + djo) — - (^j + mdjo + -mejmngmnO + -^^jmnHmri^ 



Pj 



m 
hi 



(31) 



(Note that j42^ uses a different sign convention for pj and p' than we do.) Although we 
neglected them in Section El we have included in Hpw all possible contributions arising 
from a nonzero = r° — 7*^. It is convenient to combine similar terms in (|HT|) and write 
the Hamiltonian in the form 



H 



n \ P PjPk «i R j R k 

Fw = {l- Coo) ^ - co'fc)^ + - Bjcy^ - BjkPj(y 



G 



ijk- 



m 



(32) 



The notation (jk) represents symmetrization with respect to the indices j and k, de- 
fined for example, as C(jfc) = Cjk + Ckj- Note that, in any fixed frame, Coo and the trace 
components of Cjk may be absorbed into the mass of the fermion and therefore are not 
observable; we shall henceforth ignore these terms. 

The nonrelativistic velocity operator can be computed directly from this Hamiltonian, 

as 



i [H, x] 



Explicitly, it is given by 

Vk 

where the constants a^, Bjk-, and Gjki are defined to be 



Pk Pj ak i :r P^rr^ 

m m m m 



B 



•jk 



m 



aj — m [Cqj + Cjo) — mCj, 
1 



bo^jk — {dkj + dooSjk) — mekim (^-^^mij + gmoo^ji^ — ^jkiHio 



(33) 



(34) 



(35) 
(36) 



^The same methods could be applied to the fully rclativistic Foldy-Wouthuysen Hamiltonian given in 
as well. 
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and 



G 



jki 



1 

m 



m {doj + djo) — - [bj + mdjo + -mejmngmno + -^jmnH„ 



Jkl 



mnQmnO 

) 5jk — mejlm {gmOk + gmko) 



(37) 



The result fl34|) agrees with the nonrelativistic hmit of all our other calculations. The 
same velocity may also be found by taking the direct Foldy-Wouthuysen transform of the 
relativistic velocity operator. 

Now we may determine how the modified velocity v affects the bulk motion and 
spreading of a wave packet. Suppose an initial wave packet is specified as ip{f,t = 0). 
The subsequent time dependence follows from the use of the time evolution operator: 



^(f,0) 



-iHLvtjM 



(3J 



where Hiy is the Lorentz- violating portion of the Hamiltonian and V''-^^ (r, t) is the conven- 
tional time-dependent wave packet in the absence of Lorentz violation. The decomposition 

____ 2 

(jHH|l is valid because Hlv commutes with the conventional Hamiltonian To lowest 
order in the violation parameters, the exponential e~*^^^* can be expanded to yield 



(0), 



(39) 



As an example, we consider a (p) = 0, spherically symmetric packet with its spin directed 
along the 2;-axis (s = ±1) and an initial spread determined by A. The corresponding wave 
function is 



7/'(°)(f,t) 



A 



0FfA2 + i± 

V y m 



exp 



2fA2 + z± 

\ m 



(40) 



where Xs is the two-component spinor appropriate for the spin state. 

We may calculate the action of Hlv on the wave packet using the relations 



-iVr 



7/;W(f,t) = i 



fA2 + z± 

\ m 



and 



PjPki^^^\r,t) 



1 



A2 + 



5, 



'jk 



(41) 



(42) 



This gives a probability density of 

A2 



(r,t)|2 ^ \N{t)\'exp 



{f-vHf + 



(c(ifc) - sG 



{jk)3 



m2 (A4 + 4 



(43) 
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where f| = ak/m — sBk3 is the group velocity of the wave packet and N{t) is a normahza- 
tion factor. The velocity matches the result of ()34|) for zero momentum, as expected. The 
c and G terms contribute ellipsoidal deformations in the shape of the wave packet as it 
spreads; however, in calculating this probability density, we have implicitly summed over 
the possible spin states, and therefore any spin modulations present in the wave function 
do not appear in this formula. 



4.2 Spin-Dependent Effects 

Although is s-dependent, the expression for the probability density does not tell us 
anything about the spin state of the particle. Since we have seen that there can be a com- 
plicated interplay between the velocity and the spin in Lorentz-violating fermion theories, 
it is worthwhile to look at spin effects in more detail. Our examination of wave packet 
spin structure will not be completely general; in fact, we shall consider a theory in which 
6° is the only Lorentz-violating parameter. However, this single parameter is sufficient 
to generate many interesting effects. We shall use a specialized approximation method, 
involving large time position-momentum correlations and the method of stationary phase, 
to extract information about this theory. 

The Hamiltonian governing the nonrelativistic time evolution in the 6° theory is 



2m 



-cr. 



m 



(44) 



where Up- = <? ■ p is the Pauli spinor corresponding to the direction of the momentum. 
This will act on a Gaussian wave packet, with mean momentum (p) = and the spin 
oriented along the positive z-axis. The wave function at zero time is given by 



Mr,t = 0) 




1 




1 





(45) 
(46) 



The nonrelativistic approximation requires that the length scale A satisfy mA ^ 1, and 
we shall also assume that 6° A -C 1, so that the second term in pij) may be treated as a 
small correction. 

In terms of the eigenspinors of CTp-, the wave function is 



Mr,t = Q) = (2v^A)Y 



d^p 
(2^ 




COS 1 


9 




+ sm- 


e'^ sin 1 





sm^ 
— e*"^ cos 



2 J 



(47) 
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where 9 and (j) are the spherical coordinates corresponding to the direction oip. Therefore, 
the time-evolved wave function is 



2v^A 



(27r)3 
9 
2 

d^p 



sm 



COS I 

e*"^ sin ■ 



(2vr)= 



-e^-^cosf J 



cos^ +isin^cose 
te^ sm — — sm 9 

m 



(48) 
.(49) 



The expression fl49|) is quite complicated. However, for our analysis of the velocity, 
we are primarily interested in the probability density function |'06(^, t)\'^, and that only at 
large times t ^ mA^. At such times, the different momentum modes will have separated 
themselves in space. We may therefore identify each point in space with a particular value 
of the momentum [up to a positional uncertainty of (9(A)] . However, this identification 
will not be the same for the upper (spin up) and lower (spin down) components of the 
spinor wave function 

The correct position-momentum correspondence may be found for each component by 
the method of stationary phase; writing the wave function as 



d^p 
(2^ 



^i^l {p,r,t) 
^i^2{p,f,t) 



(50) 



the position may be determined by solving Re{Vp-$j(p, r, t)} = 0. If we neglect those 
terms that are smaller than 0{t), we find 



pt 



m 



pb'-h 



™ cos2 lfi!^+sm2 M!^cos2e 







(spin up) 
(spin down) 



(51) 



If we can invert this expression, to obtain p as a function of r, we may then determine 
the spatial probability density, using that fact that the density function in momentum 
space is time-invariant and known. The dominant contribution to p{r ) is p ~ mr/t, 
and we may neglect the corrections to this expression in the argument \p\}pt/m of any 
trigonometric functions as long as t -C m/{W)'^. Making this approximation, we find 



p [r 



m / 6°t 
~ — r H 



cos 9 



t 



m cos^ h^r + sin h^r cos^ 9 



(52) 



for the spin up component. Note that p and r point in the same direction, so 9 is the 
polar angle for r as well as p. 

Our use of the stationary phase method entails three significant approximations. There 
are obviously small contributions to the wave function coming from regions of phase space 
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where $j is not stationary, and these we have neglected. We have also neglected any effects 
from the imaginary part of Vp*$j; this means, for example, that the initial uncertainty A 
cannot contribute to this part of the calculation. (A will play an important role in our 
final expression, however.) Finally, we have neglected any terms that do not grow as 0{t) 
in the large t (or equivalently, far field) limit. This means neglecting any contributions 
from Vp-^. The gradient with respect to p of 6* is a time-independent constant, and its 
coefficient will be a bounded combination of trigonometric functions. This term therefore 
does not grow linearly with t when t ^ niA/b^, and so does not represent a contribution 
to the velocity. All three of these approximations are standard elements of the method 
of stationary phase, although they manifest themselves in slightly unorthodox fashions in 
this Lorentz- violating problem. 

We also need to know the probabilities for the spin of the particle to be either up 
or down at a given time. These may be found by taking the magnitudes squared of the 
matrix factors in ()49|). Then, since each p value has a statistical weig ht of e-P^ in the 
momentum-space probability density, the position-space density is simply given by 



\iJb{r,t)\^ oc Je 



-[p(r)]2A2 



2 \p(r)\b'^t I • 2 

COS ^^^^-^ h sm 

m 

■ 2 \p(r)\b°t 
sm ^^^^-^ 



(53) 



where J = \d{pi,p2,P3)/d{x,y,z)\ is the Jacobian of the transformation from p-integra- 
tion to r-integration. We have continued to use a two-component spinor notation; the 
two matrix components of represent the probability densities for finding the particle 
at a given position with a particular spin orientation. 

The Jacobian is dominated by the r -independent term J ^ {m/t)^, but there are 
also spin-dependent corrections. These correction are obviously zero for the spin down 
component, but for the spin up portion of the probability density, they could be non- 
trivial. Sufficiently far from the origin, however, the corrections may be neglected. An 
examination of the O (b^)^ correction to the spin up Jacobian will illustrate why. Since 
all the off-diagonal terms in the spin up Jacobian matrix (the Opj/dr^ for j ^ k) are 
necessarily O {b'^)^ , and these terms cannot appear singly in any product that makes a 
nonzero contribution to the determinant, the diagonal approximation. 



J 



dpi dp2 dp3 
dx dy dz 



(54) 



is valid up to corrections that are second order in 6°. Expanding the partial derivatives, 
the leading order contribution to the spin up J becomes 



bH^ cos 9 

m cos^ b^r + siv? Wr cos^ 9 



(55) 



We may treat the trigonometric functions of b^r as constants, since their derivatives will 
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only contribute at C (6°)^ . So to first order, J is 



J 



m\3/ 26°i cos^ 
t / V mr cos^ 6°r + sin^ Wr cos^ ^ 



(56) 



The Lorentz- violating correction appears with a factor of ^cos^. At higher order, W 
will always appear either in connection with this same expression, or in (which we 
have already assumed is small). So for r » ^cos^, we may neglect the &°-dependent 
contributions to J. 

Collecting all the necessary elements, we find that our final expression for the proba- 
bility density is 



exp 



(r + 



Along the 2;-axis, the density profile is 



bPt cosg 

m cos'^ 6"r+sin^ fe'^r cos-^ 9 



(cos^ b^r + sin^ 6°r cos^ 6 



exp ^^r"^] sin^ fo°r sin^ 9 



\iljb{z,t)\'^ oc exp 



t2 



m 





' 1 ' 








(57) 



(58) 



Along this direction, the condition we have imposed on r — that f ^ ^ cos 6 — is quite 
strong. For large \z\, the exponential tails indicate the presence of an increasingly broad 
peak that moves with velocity —W/m; however, we expect that the Gaussian structure of 
the distribution will be significantly disturbed at smaller values of \z\ by the correction 
coming from J. The wave packet spreading occurs more rapidly than the 6-induced drift; 
however, it would be possible in principle (although not in practice) to measure the latter 
affect. 

In the xy-plane (where ^ ^ ^ cos 9 is automatically satisfied), the probability density 
has a significantly different form. In terms of p = ^/x^~+lp, the density is 



|^/'b(p, t)|^ oc exp 



cos^ 
sin^ Wp 



(59) 



(Note that this represents the joint probability density for x and y, conditional that z — 0; 
it is not a density for p itself.) It is in the xy-plane that we might expect to observe a 
bifurcation of the wave packet. Such bifurcation would be expected to manifest itself as 
follows: Once the distance Wt/m becomes larger than the initial positional uncertainty A, 
the peak in {ipbip, ^) P should move away from p — 0. The relevant times for this occurrence 
fall within the range of validity of our approximations, mA/b^ <^t <^ m/ (6°)^. However, 
no such effect is in fact evident. We see only unmodified wave packet spreading, along 
with a sinusoidal modulation of the spin orientation probabilities. 



17 



We see here that the 6° interaction is capable of inducing spin flips, but these flips 
only appear at O (W) 



seen in ordinary Dirac 



and higher. This situation is analogous to the Zitterbewegung 
theory. In the Dirac theory, the velocity and the spin are not 
constants of the motion; they undergo oscillations around their classical mean values. In 
the nonrelativistic, Lorentz-violating theory, a similar situation exists. If we go beyond 
leading order, then neither the spin nor the velocity are generally conserved. We have 
seen the consequences of this in our analysis. The spin flips are obviously related to the 
nonconservation of a, and the nonconservation of the velocity also has an evident effect. 

(6°)^ terms become important. 



inear function of t; this signifies 



On space and time scales large enough that certain O 
the Lorentz-violating correction to r (p) ceases to be a 
that the velocity becomes time dependent. 

Although we have seen no evidence for wave packet bifurcation in the plane perpen- 
dicular to the motion of the expectation value of the position, such bifurcation definitely 
can occur along a direction parallel to (p). An initial wave packet 



1 



C- 



(60) 



with (ps) large compared to both A ^ and 6°, will split apart under a time evolution 
governed by (jUj). The spin-up component propagates in the z-direction with velocity 
(pa)-'' 2^ while the spin-down component moves with velocity If (in contrast with 



our earlier example) WA > 1, then after a time t ~ mA/b^, the spin states will become 
well separated in space. This represents a sort of Lorentz- violation-induced Stern-Gerlach 
phenomenon. However, in the time required for the spin components to separate, the 
center of mass of the wave packet will move a large distance V^t ~ {p^)A/b^. 



5 Summary 

In this paper, we have considered the role of the velocity in Lorentz-violating quantum 
theory. In the presence of Lorentz violation, the velocity becomes a fairly complicated 
object. In order to deal with this complexity, we have used several different techniques 
in our analysis. The different methods have also been useful in our efforts to tackle some 
related subsidiary issues. Although a number of these methods, as well as some of our 
results, are fully relativistic, we have concentrated especially on the nonrelativistic regime. 

In the relativistic regime, new complications, such as potentially superluminal speeds, 
may arise. In this regard, there is a fundamental difference between the velocities asso- 
ciated with the M terms and those associated with F'^. In the presence of M only, a 
remains the velocity operator, and the Zitterbewegung-fiee portion of a has a magnitude 
that is always less than or equal to one, with equality only at m = 0. However, the F'^ 
terms do not necessarily give rise to velocities satisfying \V \ < 1, and this fact may pose 
problems for causality 3j. It seems likely that at these high speeds, some new, possibly 
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nonlocal physics will come into play, which will ensure that some form of causality is 
preserved. 

We have also seen that there are momentum and spin states that have no well-defined 
velocity, and because of this, the analysis of scattering experiments can become much 
more complicated. Cross sections for scattering between generic polarization states may 
not be well-defined. However, this does not rule out the calculation of certain particular 
cross sections. For example, in jHEl it was shown to be possible to perform all orders 
evaluations of Compton scattering cross sections, but only if a particular polarization 
basis, for which the group velocities were time- independent, was used. 

When the number of Lorentz-violating coefficients is relatively small, it may be rea- 
sonable to solve the modified Dirac equation exactly. We did this for the iP theory in 
Section A solution of the Dirac equation automatically entails a determination of the 
energy-momentum relation, from which the group velocity Vg may be determined. This 
method also involves solving for the exact propagation modes of the field, so it is easy to 
see precisely which particle and antiparticle states have well-defined velocities. 

We also used operator techniques, based on the matrix structure of the Dirac equation, 
to analyze the velocity. This allowed us to demonstrate several important facts. We 
showed that the group velocity found in Section IT^ was not merely a result of anomalous 
dispersion effects. We also demonstrated that the components of the bulk velocity V 
may be unavoidably noncommuting, so it may not be possible for the Lorentz-violating 
contribution to V to have well-defined projections along more than one axis. Finally, we 
have derived relativistic operator expressions for V to at least first order in a, 6, c, e, /, 
and H (and a similar calculation for d is located in the Appendix). 

The operator method also allowed us to address some questions about the Zitter- 
bewegung motion in these Lorentz-violating theories. We found that to leading order, 
none of the Lorentz-violating terms present in M changes the scale of the Zitterbewegung 
oscillations. 

Fundamentally, any group velocity is really a property of a wave packet, so in Sec- 
tion I3J we looked explicitly at the effects of Lorentz violation on these packets. We first 
specialized to the nonrelativistic case, introducing the leading-order Foldy-Wouthuysen 
Hamiltonian Hpy/. From Hpyy, it was a simple matter to obtain an expression for the 
nonrelativistic velocity operator and to establish its depedence on each of the Lorentz- 
violating parameters, including g. This result flMj) will be useful in all subsequent calcu- 
lations of nonrelativistic, Lorentz-violating velocities. 

We then showed explicitly how this velocity arises, by studying the time evolution of 
a particular wave packet. Not all of the Lorentz-violating coefficients contribute to the 
bulk velocity, but even those that do not may affect the structure of the packet. The 
spreading of the wave packet ceases to be spherically symmetric, because it is deformed 
by the effects of Lorentz violation. 

We also identified some interesting spin effects in our analysis. A free particle that 
begins in a definite spin state may no longer be in that state at later times. The wave 
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packet may bifurcate if the initial state is not an eigenstate of the velocity, or the spin 
orientation probabilities may develop spatial oscillations. These result indicate that there 
is a great deal of richness in the spin structure of these wave functions. 

The results we have derived should have a large number of applications in Lorentz- 
violating physics. We have provided a very general analysis of the character and role 
of the velocity in the SME, and so in any situation in which information about the 
asymptotic propagation states of this theory are required, our expressions should prove 
useful. We anticipate that this may include analyses of scattering data, decays, and 
particle oscillations. 
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Appendix: Velocity Operator in the Presence of d 

It is possible to use the relativistic operator methods described in Section |21 to study the 
velocity even when the kinetic term in the action is spin-dependent (i.e. in the presence 
of d or g). However, the fact that the spin is not conserved for a free Dirac particle makes 
these analyses less straightforward than those considered above. We shall outline here the 
calculation for the d"'^ case and display the difficulties associated with it. As in Section IHl 
we assume that d'^^ = 0. 

The presence of 75 makes the theory with d"^^ much more complicated than the similar 
c'^'^ theory. Both the single-particle Hamiltonian 

H = ajpj + dijaipjj^ + dojPjj^ + (3m, (61) 

and the velocity operator 

Xk = ak + dikaa^ + c?ofc75- (62) 

involve 75. Since they have different Lorentz structures, we shall consider dkj and d^j 
separately. 

In the presence of dkj, the equation of motion for ak is 

dfc = i {-2akH + 2pk + 2dkjPj^5) ■ (63) 
To first order in d, we may use the Lorentz-invariant expression 

75(t) = T.,p,H~^ + [75(0) - T.,p,H'^] (64) 
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for 75 in (jU^ . The equation then has the leading order solution 



akit) 



afc(O) - PkH + {dkjPj) (Tjipi) H 



-2 



-2iHt 



te 



-2iHt 



(65) 



The di^aij^ = dnJli term in Xk also contributes to 14 . The time-independent part of 
is S/(0) + ^eimnOimiS^)PnH~^ jHj. This operator represents the conserved part of 
the spin, and the Q;m(0)-dependent term vanishes in the nonrelativistic limit. So the 
Zitterbewegung-bee velocity is in this case 



Vk ^ PkH-' + {dk,p,) (I^ipi) H-^ + d, 



Ik 



{0)PnH-' 



(66) 



Because the velocity is spin-dependent, the initial conditions enter in an unavoidable 
fashion. As a result, V, like a, possesses a nontrivial matrix structure. 

We now consider a doj term. The inclusion of doj in the action produces a single- 
particle Hamiltonian very similar to the H we obtained in the presence of a W interaction. 
The only difference is that the coefficient of 75 in H is now momentum-dependent. How- 
ever, we may utilize our results from the 6° case to solve for the velocity in this situation as 
well. Let us fix k; for the remainder of this paragraph, this index is not to be summed over. 
If p is oriented in the fc-direction, and the spin is also quantized along the k-axis, then 
all our earlier operator results for 6° continue to hold. In particular, the Zitterbewegung- 
free part of is approximately PkH~' + {dojPj)T,kH-' = pkH'^il + dofcSfc), while the 
mean value of ai for / 7^ k vanishes. Since the velocity (jU^ contains an additional, 
a -independent, contribution, the total velocity along the /c-direction is 



Vk^PkH-\l + 2d^k^k 
while in an orthogonal direction / 7^ k, the velocity is 

Vi ^ doi{T.kPk)H-\ 



(67) 



(6J 



This agrees with results obtained from the energy-momentum relation, 

E"^ = 171^ + {spk + dojPjf , (69) 
where s is again the helicity; from (j^Hl) we may obtain the exact group velocities 



-^9Jk 



^9)1 



Pkjs + dgkf 

E 

sdpiPk 
E ' 



(70) 
(71) 



in agreement with (j67jl and (j68|l . 
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The difficulties with these expressions are twofold. First, there is the unattractive 
dependence on the conserved spin operator S/(0) + |e;m„am(0)p„if~^ in (jUBj) and the 
resulting matrix structure of the velocity. Second, there is the fact that (jU7j) and (jUH|) 
differ in form. We cannot combine these two equations into a single tensorial expression, 
because we broke the observer rotation symmetry with our choice of a spin quantization 
axis. While both these difficulties are surmountable, they do show that the operator 
method becomes more cumbersome in the presence of spin-dependent Lorentz-violating 
parameters. 
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